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In  this  paper,  the  steady,  irrotational ,  subsonic  flow  of  a  gas  around  a 
given  profile  is  studied  in  the  case  of  arbitrary  space  dimension  greater  than 
two.  We  prove  that  the  solution  of  this  problem  exists,  is  unique,  and 
depends  continuously  on  the  incoming  flow.  This  extends  the  previous  results 
of  Bers  and  of  Finn  and  Gilbarg. 
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SIGNIFICANCE  AND  EXPLANATION 


It  is  clear  that  the  study  of  the  gas  flow  around  a  body  moving  at 
subsonic  speed  is  of  primary  importance  in  aerodynamics*  Extensive  studies 
have  been  done  on  the  problem  of  existence  and  uniqueness  of  steady 
irrotational  subsonic  flow  of  a  perfect  gas  past  a  given  profile.  Bers  (1954) 
treats  the  plane  subsonic  flow  using  the  theory  of  quasianalytic  functions. 
Finn  and  Gilbarg  (1957)  deal  with  three  dimensional  flows  with  Mach  number 
less  than  0.7.  In  this  paper  we  prove  the  existence  and  uniqueness  of  three 
dimensional  flow  by  using  and  improving  a  priori  estimates  obtained  earlier  by 
several  authors.  Our  results  can  also  be  extended  to  higher  space 
dimensions.  Furthermore  we  allow  the  flow  to  be  arbitrarily  close  to  sonic 
speed. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 
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* 

Guang-Chang  Dong 

The  existence  and  uniqueness  of  steady  irrotational  subsonic  flow  of  a  perfect  gas 
past  a  given  profile  has  been  studied  extensively  in  the  two  dimensional  case*  In  {1], 
Bers  proved  the  existence  and  uniqueness  of  plane  subsonic  flow  around  a  given  profile. 

For  higher  space  dimensions,  few  results  have  been  previously  obtained.  Finn  and  Gilbarg 
[2]  proved  existence  and  uniqueness  in  three  dimensions  provided  the  velocity  was  not  too 
large  (the  maximum  Mach  number  less  than  0.7). 

In  this  paper  we  use  the  idea  of  [2]  together  with  an  improved  a  priori  estimate 
(extending  the  method  of  (4) )  to  prove  the  existence  and  uniqueness  of  the  solution  in  the 
three  dimensional  case.  We  also  extend  the  result  to  higher  dimensions. 

In  the  following  we  always  suppose  the  dimension  of  space  is  n(n  >  "M  ,  and  use  the 
summation  convention 


n 

ax.  =  T  ax  , 
ii  ,  i  i 

i=1 

and  denote  the  vector  (x1,...,xn)  by  x. 

The  steady  irrotational  gas  flow  in  n  dimensional  space  can  be  described  by  the 
/elocity  potential  vA(x)  ,  satisfying  the  equation 


S2-  (p  -  o  <i> 

o  X .  v  oX 
1  1 

where  p  represents  the  density  of  gas,  which  is  a  given  positive  function  of  velocity 
q,  where 

/TT  3v 

q  =  / lu  .  ,  u  .  =  — 

l  l  Ox  . 


( 1 )  can  be  written  as 


3  y  _ 

ai j  Ox . Ox . 
i  0 


U'  ' 
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(3) 


Because  (1)  and  (2)  are  rotational  invariant,  the  eigenvalues  X . ,  X  of  the  quadratic 

i  n 

form 

aijai°j 

can  be  obtained  by  letting  u1  =  q,  u2  =  •••  =  un  »  0,  i.e. 

X,  »  p  +  p'q,  X,  =...*  X  =  p 
i  2  n 

So  that  if 

p  +  p'q  >  0  (4) 

at  this  point  the  quasilinear  equation  is  of  elliptic  type,  and  the  flow  (1)  is  subsonic  at 

this  point.  If  in  some  region  S2  we  have 

inf ( p  +  p'q)  >  0  , 

S2 

then  the  flow  (1)  is  a  subsonic  flow  in  fi. 

Assume  that  the  function  p  =  p(q)  is  a  positive,  three  times  continuous 
differentiable  functions,  defined  in  the  interval 

°  <  q  r  qUm  4  -  , 


and  assume  that  a  number  q  exists  (we  called  it  the  critical  velocity)  0  <  q  <  q 


such  that 


p  +  p'a  =  — -  (pc!  >  0  (0  <  q  <  q  ) 

a  a  c 


lim' 


(5) 


For  definiteness,  assume  that  when  q  >  oc ,  the  left  hand  side  of  (4)  is  non-positive. 
Hence  (2)  is  subsonic  if  and  only  if  0  <  q  <  q^.  We  also  assume  that 

p'(0>  =  0  .  (6) 

In  gas  dynamics,  the  density  of  a  perfect  gas  is  given  by  the  dimensionless  formula 


o  =  b  -  <i 


,  V-1 


<1  <  y  <  2,  y  is  a  constant) 


o'b  “  ' 1  ■  q2)  b  ■ 


IZY 

Y-1 


Y+1  2, 


hence 


qc  “/y+r  qfim 


y-1 


The  flow  is  culled  subsonic  when  0  <  q  <  q  ,  which  agrees  with  the  discussion  above. 

c 

The  profile  T  is  a  bounded  closed  surface  in  n  dimensional  space  (it  can  be  some 
closed  surface  also) .  Assume  that  a  constant  exists  such  that  0  <  <  1  and 


(2+V 

r  e  c  (7) 

The  region  outside  T  (not  containing  the  point  at  infinity)  is  denoted  by  ft.  We  only 
consider  those  satisfying  the  condition:  any  closed  curve  in  ft  can  be  deformed  into  a 
point  without  touching  T . 

The  flow  around  T  means  tnat  the  solution  of  (1 )  in  ft  satisfies  two  conditions  as 
follows.  The  boundary  condition 


3£i 

3N  | 


=0  (N  is  the  interior  normal  of  ft  on  D 


and  the  condition  of  uniform  flow  at  ft  (called  uniform  incoming  flow) 


grad  f  I 


=  u  (u  is  a  constant  vector) 


(8) 


(9) 


Usually  for  the  problem  of  flow  around  a  given  profile  we  give  the  magnitude  and 
direction  of  incoming  flow.  Without  loss  of  generality  we  can  assume  that  the  direction  of 
incoming  frow  is  the  positive  x1  axis  (otherwise  after  an  axis  rotation  we  can  reach  this 
situation),  i.e. 


00 

u  =  (U,  0,...,0)  the  const  U  >  0  is  given  •  (10) 

A  slightly  different  problem  is:  Given  the  direction  of  incoming  flow  but  not  the 
magnitude,  i.e. 

u  *  (U,  0,...,0)  U  >  0  is  not  given  (11) 

and  given  the  maximum  value  of  speed  instead,  i.e. 

sup  q  *  Q  is  given  •  (12) 

a 

We  shall  study  the  existence  and  uniqueness,  find  the  flow  and  other  properties  of  the 
solution  of  the  above  two  kinds  of  flow  problem  in  the  subsonic  range,  i.e.  when 

Q  <  qc  •  (13) 

§1.  Preliminary  study  of  the  linear  problem 
Consider  the  following  problem 


-  3- 


2 

L<y)  -  Vx)  aHr +  Vx>  fr  ■ f(x)  (x  e  n'  bij  -  v 


(i.D 


3f 


3n(  ”  *o  *N  iS  tbS  ^nter*or  normal  of  ft  on  T) 


(1.2) 


*1  *  0 


(1.3) 


where  the  given  functions  b^,  bA,  f  c  C(0,(ft>,  fQ  c  C ( 0 1  ( T ) .  and  a  positive  constant  a 


exists  such  that 


-  2  .  .  ,1.2 
ola,  4  b.  n  a  <  —  La  . 
i  1)  i  ]  o  1 


(1.4) 


Further  conditions  on  b^ ,  b^,  f  are  given  below.  Assuming  that  the  solution  of  problems 
(1.1),  (1.2),  (1.3)  exists,  and  we  deduce  some  a-priori  estimates. 

Assume  that  a  fixed  small  positive  number  T  such  that  t  4  tq  (the  definition  of 
tq  see  (7))  exists  such  that 


lb  (x)  -  b I  <  K  r  (x  e  ft) 

13  ^  J  IX 

|bij(x)  "  bij(V,l  *  K2rxy|x“y|T  (X'  V  E  n  5 

|b. |  4  K ,r"1*T 

1  3  x 


(1.5) 

(1.6) 

(1.7) 


where  the  constants  K ^ ,  Kj,  ...  depend  on  o,  t,  T  only.  Moreover,  assume  f  is  a 
bounded  function  satisfying 


I  f |  4  K  max  | f | r  n  T  (1.8) 

4  ft 


where 


r  »  d(x,  D,  r  =  min(r  ,  r  ) 
x  xy  x  y 

d(x,  P)  is  the  distance  between  x  and  T.  We  have 
Lemma  1 . 


(1.9) 


Proof.  Let 


If |  4  K  (max  If  |  +  mew  |f|] 
T  0  ft 


f  -  max  |f  |  +  max  |  f  | 
T  ft 


(1.10) 


(1.11  ) 


-4- 


The  uniqueness  of  the  solution  of  (1.1),  (1.2),  (1.3)  is  easily  deduced  from  the 
maximum  principle  [5].  Thus  for  =*  0  we  have  ¥  =  0,  and  in  this  case  (1.10)  is  true. 


When  H*  ¥  0,  let 


L-(*>  ■  V”’ 


R  ■■  /B  ..(•»)  x  .  x  . 

i  j  13 


B  («>)  ig  the  algebraic  complement  in  the  matrix  (b.  ,(“>)).  We  have 
ij  13 

2-n-  j  -n-  | 

lJr  ]  -  j  ( n-2+  ~)R  .  {  3.14) 

When  R  >  K  ,  where  K,  is  a  sufficiently  larqe  constant,  from  (1.5),  (1.7),  (1.8), 

O  ® 

(1.14)  we  have 


2-n-  -  -n-  — 

L(R  2)  -  ^  =  “  <n-2+  j)R  2  +  0 ( R~n_T )  >  0  . 


Without  loss  of  generality  we  can  assume  that  the  origin  lies  inside  r,  then 

min  R  >  0.  Take  a  =  a(cr,T,n  comparatively  large,  by  using  (1.5),  (1.7),  (1.13)  we  have 
xcii 

UR*01)  »  a(a+2)R'“~4b.  .B.  .  (»)x,  B.w(«.)x  -  aR~0,'2b.  .  B .  .  ( »)  -  (1.16) 

i]  ij  k  ih  h  11  ii 


<*'Cr2biBik(-,1,k  >  K1R""'2 


From  (1.15),  (1.16)  we  can  select  a  constant  such  that 


2-n-  ~ 

L(R  +  KgR  "‘j  -  —  >  0  ( x  e  71) 


1  n  2'n_  n 

T  |  =  V  +  f°(R  +  KflR  a)  . 


From  (1.1),  (1.17),  (l.lfl)  we  have 


L(4»  )  <  0,  L('i»2)  >  0 


-5- 
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From  (1.3)  we  have  ¥  (®)  **  0,  hence  min  v  <  0,  this  minimum  can  only  be  taken  at 

a 

»  or  on  T  by  (1.19).  Similarly  max  ¥  >  0,  and  this  maximum  can  be  taken  at  »  or 

ft  2 

on  T*  Hence  we  have  four  cases  as  follows. 

Case  1.  min  ^  =  ^(Pg)  <  0,  PQ  e  T;  max  ^  -  '*'2(Q0)  >  °'  e  F  * 

( 2 ) 

Extend  the  method  of  (6)  to  prove  that  (1.10)  is  true.  Because  T  e  C  ,  there 

exists  a  positive  constant  y  ®  y(F),  such  that  for  any  point  P  e  T  we  can  draw  an 

exterior  tangent  sphere  with  radius  y  which  lies  entirely  inside  ft  except  point  P. 

Draw  the  exterior  tangent  sphere  with  radius  y  at  Pq,  denote  its  center  by  P^. 

Consider  the  function  _ 

-hPP2  2 

¥  (P)  *  Y^P)  -  min  V  -  kV  (e  -e  Y  )  (P  e  ft)  (1.20) 

ft 


where  k,  h  are  positive  constants  we  may  choose  at  our  disposal.  Take  h  »  h(o)  large 
enough,  from  (1.17),  (1.20)  we  have 

LOfj)  <0,  |  <  P^P  <  y  .  (1.21) 

From  (1.20)  we  have 


If 


Yj  >  0,  when  P^P  =  y  and  ( PQ )  =0 


(1.22) 


Y3  >  0  when  P^  =  ^  (1.23) 
is  true,  then  from  (1.21),  (1.22),  (1.23)  we  see  that  min  Y^  <  0  when  ~  <  P^P  <  y  can 
not  he  true,  hencp  min  =  Y^tP^)  =  0.  From  (1.2),  (1.1B),  (1.20)  we  have 


iy 

3 

3N 


Ip  Vpo' 

0 


!s<r 


2-n- 


v‘d,| 


-  2YkhY0e-hY 


>  0 


(1.24) 


Pick  the  constant  kQ  so  that 

,2  -  2-n-  \ 

2yK<,he~  Y  >  1  +  max|-^(R  *  KRR~“)  | 


1 


then  (1.24)  can  not  be  true  when  we  take  K  *  K, 


hence  (1.23)  can  not  be  true,  in  other 


words,  we  have  P2  satisfies  p1p2  “  2  such  that  <  0*  From  (1*20)  we  have 


S’  1  ( P2  )  -  min  <  K^Y0,  i  <  d(P2,  PQ)  <  I1  .  (1.25) 

ft 

Similarly  there  exists  a  Q2  satisfying 

V  1v 

(1.26) 

Because  of  (1.6)  we  can  apply  Harnack's  Inequality  [7]  to  the  nonnegative  function 


max  f2  -  ^2 (Q2 )  <  K11f°/  j  <  dffi2'C0)  *  f1  * 


-  min  Y 1  and  obtain 


'f  < C_ )  -  mm  Y  <  K  [Y  (P  )  -  min  Y  +  max|L(Y  )|] 

J  s  *  a  n 


Combine  (1.11),  (1.18),  (1.25),  (1.27)  we  have 


f1(C2)  -  mi"  «  K^f 


(1.27) 


(1.28) 


From  (1.18),  (1.26),  (1.28)  we  have 


ma  i*  -  min  ¥  <  max  4*  -  min  4*  <  (K 

ft  ft  ft  ft 


K13>* 


From  (1.3)  we  have  max  Y  >  0,  min  Y  <  0,  so  combining  the  above  expression  and  (1.11)  we 

n  (5 

obtain  (1.10). 


Case  2. 


min  Y 

(5 


Y  («)  =0,  max  Y 

Q 


VV  >  o,  q0  t  r 


By  lemma  2,  there  exists  positive  constants  R  (o,T,D  and  K14,  such  that  the 

Harnack  Inequality  for  the  positive  function  R2  "v,  (the  definition  of  R  see  (1.13)) 


for  any  R1  >  RQ  holds.- 


msx  (R2  "f,  )  <  K  min  ( R2  nY  )  +  K  max  |  f  | 
R=Rq  1  14  R=R,  1  15  0 


(1.29) 


Proof  of  (1*29)  see  lemma  2.  From  (1.29)  we  have 

max  Y  =  r"'2  Max  (R2'nY))  <  K^rJ'2  max  (R2'^)  +  K^R^'2'*'0 
R=R0  R=R0  R=R, 


R  n-2 

( 

'M'-R 


— )  max  Y. 


1 


K  r"-V 

K15R0  Y 


(1.30) 


From  (1.19)  and  4* 


=  V 


1 


=  0, 


combining  u  .th  (1.18)  we  have 


max  f.  *  max  f  <  ma  .  4*  <  max  4*  +  K  Y 

*  2  2  I  To 

R=Ri  R=Ri  R=Ro  R=Ro 


(1.31) 


so  combining  (1.30),  (1.31)  we  have 
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R„  n-2 

r  O' 


R„  n-2 

0  i 


11  -  *14^  >  h  «  '*14^  *16  +  V  V 


,K0\  1 

Take  R-  such  that  K  I — I  =  *r  we  have 
1  14  R.  2 


max  f  1  <  K  l?  .  M . 

R=Ro 

(1.32)  is  similar  to  (1.25).  We  can  get  (1.26),  (1.27),  (1.28),  (1.10)  by  the  similar 
process  as  in  case  1. 

Case  3.  min  »  =  *  (P  )  <  0(P  e  D ,  max  f  =  t  («)  =  0. 

£2  3  2 

Obtaining  (1.10)  is  similar  to  case  2. 

Case  4.  mm  'f  =  0  =  !'  (“),  max  ¥  =  0  =  '#(«•). 

8  1  1  3  2  2 

Prom  (1.18)  we  have  ¥  =0,  contrary  to  the  hypothesis. 

Lemma  1  is  thus  proved. 


Lemma  2.  Let 


*,  >  0  <*  e  n>  • 


Prove  that  (1.29)  is  true. 


Proof.  Assume  that  the  region  R  >  R  (the  def  init  ior.  of  R  s*  e  (1.13))  lies  m 
Without  loss  of  generality  we  can  assume 

max J  f |  =  1  (1.34 

5 

b.  .(«*)*  <5.  .  .  (1.35 

13  13 

Otherwise  after  a  linear  change  of  independent  variables  and  multiplying  the  unknown 
function  by  a  suitable  constant  we  get  (1.34)  and  (1.35).  Hence  R  =  |x|.  Apply  the 
inversion  transformation 

x  . 

X.  =  — —  (i  =  1,2,...,n)  (1.36 


.  n2_n'  2v 


(1.1)  becomes 


L(v)  =  B.  ,vv  Y  +  BuVv  +  Bv  =  F  (  ( X  |  < 

r3  XiXj  i  Xj.  -13 

Notice  that  the  matrix  B^  is  not  the  same  as  in  (1.23).  From  (1.4)  and  (1.36)  we  see 


-8- 


I 


iii.  There  exists  constants  K21,  K^2  such  that 

K 

K+(X.  Y)  >  K21.  KJX.  Y)  4  K22  (IXI  «-§*,  ivl  =  K20)  .  ,1.49) 

If  K+(X,  Y ) ,  ( X ,  Y)  exist ,  then  combine  (1.44),  (1.45),  (1.47),  (1.48),  (1.49), 

and  apply  the  method  of  15],  (c.f.  in  $4  of  [5])  which  extends  theorem  1  to  the 
inhomogeneous  equations.  This  is  the  result  that  we  need.  Hence,  positive  constants  ^23* 
K24'  K25  ex*st  such  that 

K 

K,,v<0)  -  K  <  v(X)  <  K'l  v(0)  +  K  (  |  X  |  <  — §2-)  ,  (1.50) 

23  24  23  25  3 

(1.29)  follows  from  (1.34),  (1.37),  (1.50). 

We  can  construct  K+  and  K_  by  only  altering  the  method  of  [5)  a  little  as  follows. 


H  =  (K2  -X2  )  fB*  .  (y)  (X.  -  Y\)(X.  -  Y .  )  )  2 

20  1]  1  i  3  3 

where  B*  .  are  the  algebraic  complements  of  B^.  Take  K_  =*  f_(H),  then 

L(K  )  =  B.  ( X ) H  H  ,f"  +  Af 1 )  +  Bf 
11  X .  X  - 
l  1 


A  =  f  B.  .  (  X  )H 


ij(X,HXiXj  +  BiHXil/(Bij(X)HXiHx 


Apply  the  estimate  method  of  [5]  (the  method  of  obtaining  (30)  in  it)  and  (1.6), 
(1.39)  and 

f  -  1  I  '  1 

B.  =  0  (  |X|2  )  =  0  ,  |  X-Y  |  -  X20l  2  ) 

deduced  from  (1.40),  we  have 

-  -  1 

I  A  |  <  K2(.<|X-Y|n~1  +  T  +  I  X-Y  |  n  I  X-Y  |  -  K2q|2  )  . 

And  then  determining  f  etc.  (these  steps  are  similar  to  [5]),  at  last  we  get  K_(X,  Y) 
satisfying  (1.47),  (1.4R),  (1.40). 

Through  the  t ransforma t ion  K+  =  (1+2  |X|2)K*,  (1.47)  becomes 

L* ( K*  )  >0 


Applying  ( 1.45)  we  have 


L*  ( v )  =  R .  .  v  +  B#Vx  +  B*v 

11  i  j  1  i 


B*  *  1/  1  +  2|X|2)/(1  +  21x1  2 )  >  0  . 
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Let  K+  »  f  +  ( H )  •  the  selection  of  f+  Is  similar  to  [5].  The  proof  of  the  existence 
of  K *  such  that  K*  satisfying  (1.47),  (1.48),  (1.49)  is  similar  to  K_.  Hence  the 
proof  of  lemma  2  is  finished. 

Lemma  3 .  We  have 

Ifl  <  K  mind,  r2  n  Umax  If  |  +  max|f|)  (1.51) 

2  r  °  5 

under  the  conditions  (1.4),  (1.5),  (1.6),  (1.7),  (1.8). 

proof.  Because  of  lemma  1,  we  only  need  to  prove  the  validity  of  (1.51)  when  r 


large  enough.  From  (1.13),  (1.14)  when  R  >  (*28  large  enough)  we  have 


L( R2_n  -  R 


,  T 

2_n'  2( _ f_ 

*5f 


T 

-n-  — 

-  -  -  j(n  -  2  +  j)R  +0  (R_n_T)  <  0  . 


Hence  the  function 

2-n-  ^  2-n-  j 

*5(R2_n  -  R  ,/(K28n  "  K28  ’’t'  1  (1.52) 

can  not  take  on  a  neqative  minimum  value  in  Kjg  <  R  <  «».  From  (1.3),  (1.10)  we  see  that 
when  R  -  «>  and  R  -  K2g  (1.52)  is  nonnegative,  hence  (1.52)  is  nonnegative  in 
<  R  <  <».  This  proves  the  lemma. 

I k !  (  k) 

Denote  the  k'th  order  derivatives  of  F  by  D  F.  When  Fee  ((?)  or 
C(1<  +  T,(n),  let 


M  (F>  =  lub  rm+k|D(k)f(x)  | 
m,k  „  x 

xeft 


(1.53) 


M  IF)  =  lub  r”  |0(  ’f(x)  -  Dl  F(y) l/lx-yl T 

m,k+T  _  xy 

x.veft 


(1.54) 


where  m  is  a  nonnegative  integral.  For  the  meaning  of  symbols  rx,  rxv  see  (1.9).  The 
symbol  of  lub  in  (1.53),  (1.54)  means  first  to  take  upper  bound  for  all  k’s  order 
derivatives.  Let 


I  .  =  V  M  .(f),  lf|  ,  =  Ifl  .  +  M  (f) 

m,k  m,k  m,k  +  x  m,k  m,k  +  x 

i=o 
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When  H<(x)  e  C(k)(0)  and  I?  I 


is  finite,  denote  by  f(x)  e  C 


Similarly  we  define 


Lemma  4.  Assume  that  (1.4),  (1.6),  (1.7)  and 


Ib^x)  -  b.(y>|  <  K29|x-y|  Vy"T 


are  valid,  and  assume  that 


Let  ¥  e  C^+TN^)  he  the  solution  of  (1.1),  (1.2),  (1.3)  and  ^  c  C  _  ,  then  we  have 

n-  2 , 0 

f  £  Cn-2 ,2+T '  and 

|y|n-2,2+T  «  K30[Mn-2,0(4,)  +  ^'n.T1  *  <1‘57' 

Proof.  Prom  the  interior  estimate  [9]  we  obtain  (1.57). 

-  -  I],) 

Consider  a  bounded  region  satisfying  ft  o  ftQ  ^  T.  when  ¥  c  C  (ftQ)  or 

c  C(k+T) <flQ)  >  let 

M  (*>  “  lub  |D(k’!i'(x)  j,  M  It)  -  lub  |D(k),Kx)  -  Dll!)1(y)  )/|x-y|T  , 

xenQ  x,yE£l0 

=  jo  Y”'  '’’m  “  IVIk  +  V?>  * 

When  Y(x)  E  Ctk'(SlQ)  and  |Y|^  is  finite,  denote  by  fix)  e  C®,  similarly  we  define 


Concernng  the  function  defined  on  T,  let 

»»-».  -  lub  lub|D<h)(f  )|,  l¥  I  -  IV.  +  lub  lD(k),fnlx) 

0  *  u^i  _  r  0  0  k+T  Ok  r  0 

h«k  xe"  x,yeT 

'  D(k)4'  (y)|/|x-y|T 

k  ( k )  u 

where  in  lub  the  D  ^(x)  and  D  ^(y)  restrict  to  the  derivatives  of  the  same 

pa rameter. 

Lemma  5.  Apart  from  the  hypothesis  of  lemma  4,  let  b^,  b^,  f,  yQ  satisfy  the  following 
conditions 


lb.  .(x)  -  bi,(y)|  <  K31|x-y|T  (x,  y  e  ilQ) 
|bi(x)  -  bi(y)|  <  K32lx~y|T 


TQ  e  C(1  +  T)  (D  . 

And  let  T  c  C(2+t)(Q)  be  the  solution  of  (1.1),  (1.2),  (1.3),  then  we  have  T  e  C° 


and 


Proof.  See  [10] 
Let 


|T‘2+,  4  K33fM0<y)  +  1  f  1  t  +  "oW 


M  (T)  =  lubmaxd,  rm"fV: )  I  D(>< !  T(x)  | 
m,  k  _  x 

xen 


M  ,  =*  lub  max(  1 ,  r  )|D  T(x)  -  D  'ft  y )  I  / 1  x-y  I 

m.'c+T  „  x 

x.yeO 


(1.61) 

2+t' 

(1.62) 


iti  ,  =  T  m  .(T),  it i 

m,k  ^  m,j  m,k+T 


ITI  .  +  M  ( T) 
m,k  m,k+r 


Summing  lemmas  3,  4,  5  up,  we  have 

Lemma  6.  Under  the  conditions  (1.4),  (1.5),  (1.6),  (1.7),  (1.8),  (1.55),  (1.56),  (1.58), 

(1.59),  (1.60),  (1.61)  let  T  e  C(‘b(2)  be  the  solution  of  (1.1),  (1.2),  (1.3),  we  have 

IT  I  „  ,  <K,.(Ifl  +  IT*  II .  )  .  (1.63) 

n-2,2  +  t  34  n,T  0  1  +  t 

Theorem  1.  If  the  conditions  of  lemma  6  are  satisfied  then  the  solution 
T  t  C(0,(?i,  n  C(2,(a,  of  (1.1),  (1.2),  (1.3)  exists  and  is  unique.  T  e  C<2+T>(f!),  and 
satisfies  the  estimate  (1.63). 

Proof.  The  uniqueness  of  the  solution  follows  directly  from  maximum  princiole.  The 
existence  of  the  solution  can  be  obtained  by  applying  the  continuity  method,  i.e.  to  solve 
the  equation  with  parameter 

((1-9)L  +  9L)T  =  9f  (0  <  0  <  1)  (1.64) 

00 

with  boundary  conditions  (1*2),  (1.3),  where  the  definition  of  operator  L^  is  defined  by 
(1.12).  When  9=0,  through  a  linear  transformation,  the  equation  (1.64),  (1.2),  (1.3) 
becomes  the  Laplace  equation  in  an  exterior  domain  with  oblique  derivative  aiven.  The 
existence  of  a  solution  f  and  ^g  C^+T^(fi)  can  be  obtained  bv  the  method  of  integral 
equations,  see  [S],  Applying  lemma  6  to  extend  the  solution  hv  increasing  the  parameter 
b  unt  i  1  b  ;=  1 . 

The  theorem  is  thus  oroved. 
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1 


Hence  concerning  the  proof  of  theorem  2,  without  loss  of  generality  we  can  assume  the 
quasilinear  equation  (2)  is  uniformly  elliptic  when  0  <  q  <  »,  i.e.  there  has  positive 
constant  u  •  o(Q)  such  that 


2  12 
oEa.  <  a  (Di)a.  a.  <  —  Ea. 
i  ij  i  l  o  i 


Let  the  function  fix)  satisfy 


e  c(1+T)<n+n 


D*(«)  -  (0,  0, 


,0) 

'  0 


(2.4) 


(2.5) 


(2.6) 


|  D  x  >  -  D^j(®)  I  <  c^r 


-T 

X 


I D^( x)  -  D^(y)|  <  C  mind,  r  T)|x-y|T 
*  xy 


where 

T  “  T ( Q , r )  <  T  (the  definition  of  see  (7))  (2.7) 

is  a  positive  constant  which  shall  be  determined  later,  U  is  a  positive  constant  not  to 
be  fixed,  (for  the  meaning  of  rx,  rxy  see  (1.9)),  Cj,  C2  are  positive  constants.  The 
norm  of  tf(x)  is  defined  by 

lv?l  “  |U|  +  maxl^-UxI  +  inf  C  +  inf  C  . 

a 

It  is  easy  to  see  all  functions  c  form  a  Banach  space  E. 

We  wish  to  find  the  solution  of  the  following  equations: 

V0**  -  0 


(2.8) 


i  1 


15|r  -  •cos(N'  Vr 


fl  =  o  . 


Let 


*(x)  =  6QlT(x)  +  x  ]  /max  |  grad[f(  x)  ♦  |  (0*641) 

a 

where  the  definition  of  a  see  (12). 

By  theorem  1  we  have  4  c  E,  hence  a  functional  in  E  is  defined  as  follows 


(2.9) 


(2.10> 


(2.11 ) 


4  *  T(-r,  0)  (0*9*1) 


(2.12) 


I 


Now  use  Leray-Schauder  degree  theory  (11]  to  prove  that  for  any  0  <  9  <  1  the 


solution  of  functional  eouation 


*>- T(V!  ,  9)  -  0 


exists.  By  theorem  1  and  (2.8),  (2.9),  (2.10),  (2.11)  we  see  that  when  9  =  1  the 
solution  of  (2.12)  is  a  solution  of  (1),  (8),  (9),  (11),  (12),  (13). 

Let  us  verify  the  three  conditions  of  Leray-Schauder  degree  theory, 
i.  From  theorem  1  we  have 

1*1  •>  ^  •  (2.14) 

n-2,2+T  34 

Hence  T  is  completely  continuous.  T  is  also  closed,  i.e.  when  ^  e  E,  >f  e  E  and 

I  -  'ft  ♦  0  (k  •»  <•)  the  corresponding  solutions  satisfying  14^  -41+0.  If  this  is 

not  true,  i.e.  there  is  a  subsequence  {k* }  and  positive  constant  S  such  that 

14.  ,  -  41  >  6  .  (2.15) 

k* 

When  k'  is  large,  from  (2.14)  we  have 


I*.  ,  I  „  ,  <  K  (1  +  1^1)  , 

k1  n-2,2+t  34 


hence  {*  , }  is  compact  in  1...I 
{k" }  of  {k1}  such  that 


n-2,2  +  j 


i.e.  there  exist  and  a  subsequence 


«v  -  **« 


From  a  (DV) 

1  3 


Similarly  we  have 


n-2,2  +  ^ 


0  taking  limits  we  have 


aij<D^  °  • 

i  j 


—  |r  «  -cos(N,  Xl)lr,  *.| 


Hence  from  the  unioueness  of  the  solution  of  (2.R),  (2.9),  (2.10)  we  have  -  Y,  hence 

«$,  N  -  ♦  0  . 

k" 

This  contradicts  to  (2.1^),  therefore  T  is  closed. 

From  (2.11)  we  see  that  T  is  continuous  uniformly  in  9,  hence  the  degree  of  (2.13) 


is  indenendent  of 

li.  When  R  *  0 ,  <t>  *  0 ,  i.e.  (2.13)  has  only  the  solution  =  0.  Therefore,  the 


legree  of  (2.11)  is  egua1  to  1. 


I 


iii.  To  prove  the  solution  of  (2.13)  is  bounded  in  E,  in  other  words,  to  prove  that 

the  solution  of  (1),  (8),  (9),  (11)  and 

max  q  -  0Q  (0  <  6  <  1 )  (2.16) 

£2 

is  bounded  in  E,  which  shall  be  proved  in  the  following  two  sections,  first  notice  that 
from  (2.16)  we  have 

max  q  <  Q  (2.17) 

£2 


§3.  The  interior  and  boundary  estimates  of  solution. 

Let  ^  £  E  satisfy  (1),  (8),  (9),  (11),  (2.17).  Consider  (2)  as  a  linear  equation  in 
•f,  from  (2.5)  and  the  Schauder  estimate  (notice  that  we  have  (7)  and  (2.7))  we  have 

VJ  £  C(3+T>  (£2)  n  C(2+T)  (£2+r )  . 


Differentiate  (1)  about  x.  we  have 


3^ 


3  3  u 

r —  (a.  t — )  =0  (u  is  the  abbreviation  of  - — 
dx .  13  3x.  3x^ 

1  1  h 


V 


where  the  definition  of  a^  see  (3). 

Assume  the  sphere  I  x-xQ  (  <  p  lies  in  52  and  0  <  6  <  1.  Let 

1  |x-xQ|  £  p(1-0) 

C(x)  =  v(x,  p,  6)  =  \  ^  (u  -  |x-x0l)  p ( 1-0  )  <  I  x-xQ  |  <  p 

0  I  x-x  |  >  p  . 

2 

Multiply  (3.1)  by  t,  (x)u  and  integrate,  we  have 

|<;2a.  .u  u  dx  =  -2ja.  ,uu  tt  dx 

1]  X  ,  X  .  ‘  1]  X  .  X  . 

i  3  J  1  ] 

From  this  formula  and  (2.4)  we  have 

o/c2(grad  u)2dx  <  L(/ (EC2 ( grad  u)2  +  ^  u2(grad  c)2]dx 


this  and  the  following  sections  the  constants  L ^ ,  L2 , 


(3.1) 


(3.2) 


depend  on  Q,  T  only. 


an  arbitrary  positive  constant.  Take  eL^  =  —  we  have 

/ i;2(grad  u)2dx  <  L^/u2(grad  c)2dx 
Take  0  =  ^  from  (3.2),  (3.3),  (2.17)  we  have 

/  (grad  u)2dx  <  L  pn  2 

l*-*0,*2 


(3.3) 


(3.4) 
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Let  n  be  a  constant.  Multiply  (3.1)  by  the  function 

.2  . 


{ 


c  ( x)  (u-n>  (u  >  n) 
o  (u  <  n) 


and  integrate,  we  obtain  similarly 

/  t^tgrad  u)2dx  4  L /  (u-n)2(grad  C)2dx 
u>n  u>n 

We  denote  the  sets  u(x)  ?  n  in  sphere  Ix-x^l  4  \i  by  A  ,  B  respectively 

^  n  t  u  n  §  m 

From  (3*2)/  (3*5)  we  have 

5 


/  (grad  u)^dx  <  L  mes  A  — r  max  [u(x)  -  nj^ 

A  5  n,u  (6  j2  xeA 

n,u-«u  n,u 


(3.5) 


(3.6) 


Similarly  we  have 

f 


B 


2  1 

( grad  u)  dx  4  L,  mes  B  - 

6  n.u 


max  (h-u(x)] 


(3.7) 


(6p)  XEB 


n,u-0u  "  n,u 

From  (3.4),  (3.6),  (3.7)  and  applying  the  result  of  (3):  There  exist  constants  L,  and 


r  =  T  < Q , T )  (0  <  T  <  1)  such  that  for  any  x,  y  e  SI  we  have 

_T1  T1 

|u(x)  -  u(y)  |  4  L  r  |x-y| 

/  xy 


(3.8) 


Note.  In  (3)  the  inequality  (3.8)  is  proved  under  the  restriction  |x-y|  4  L  r  .  When 

o  xy 

| x— y (  >  Larxy,  the  validity  of  (3.8)  is  true  by  (2.17). 

We  now  turn  to  the  boundary  estimate.  From  (7)  we  see  that  there  is  a  local 
parametric  expression  for  a  little  part  of  T 


*i  =  VS . S-11  E  C 


(2+V 


in  the  neighborhood  of  any  point  e  T . 

Choose  variables  (£„,•••»£  )  suitably  such  that  inside  the  small  sphere  u>  with 
1  ,  _  n 

<2+V 

center  PQ  we  have  5  e  C  .on  point  PQ  *S '  *  ’  *  'S'  for  a  unit  orthogonal  axis 

system,  and  £  coincides  with  the  direction  of  normal  on  w  n  P.  The  method  for  choosing 

is  as  follows:  First  take  £.,...,£  .  such  that  they  form  the  unit  orthogonal  axis 

i  n- 1 

system  on  PQ.  Penote 


cos(N,  x  )  *  A .  (£  , 
i  l  1 


,Z  J  ( i  = 

n-1 


on  F.  Then  the  inverse  function  £ 


of 
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*1  *  V51 . C„-1l  +  en  i  VC1  +  °15n . 5n-1  +  Vl5n)dV 


satisfies  all  our  requirements. 


From  (3.9)  we  have 


=  !Li!^ 

Jjk  3x^  3x^ 


ajkir  -  Sk  <>- 11  - 1 . n> 

0 


x.  I  *0 

3n,unr 


2  r  2  -  -  3v»  , 

q  =  Eu.  »  a  u  u  (u.  «  -rt— ) 
i  ]k  3  k  3 


Equation  (1)  is  the  Euler  equation  of  the  variation  problem  6  /  F(q)dx  ■  0  where 
F(q)  =  /  pqdq.  Under  the  transformation  we  have 

. 3( x  . . . ,x  )  , 


IQ  \  A.  •  •  •  ,  A  /  . 


Tt—  (Jpa.  ,u. )  *  0 
H..  13  3 

Differentiate  (3.12)  with  respect  to  E  we  have 

n 


3  3uh 

■jj-(JB.s  +  B.  )  =  0 

35j  13  1 


B .  .  *=  pa,  ,  +  —  a . ,  a  .  u,u 
ij  lj  q  lk  k  m 


3  *  km  —  —  — 

8,  =  T7—  (Ja.  .  )pu.  +  2Jp'qa.  .  -rr —  u  u  u 

1  3£h  13  3  13  9ih  3  k  m 


From  (3.10),  (3.14)  we  have 


S  .  ,  l  =  (pi  .  ,  +  —  u  u  .) 

131P0  13  q  1  3  PQ 

Hence  when  the  radius  of  to  is  sufficiently  small,  form  (2.4),  (3*16)  we  have 


(3.17) 


I  «  hWi « !  Ex?  • 

Denote  the  image  of  u,  u  n  (!,  u  n  T  by  w,  8,  F.  (8)  becomes 

u  l  -  0  .  (3.1 

n'F 

Let  5(C<  P,  0)  be  the  function  defined  by  (3.3)  with  x,  x^  changing  to  5, 
restrict  in  this  case  that  the  sphere  with  center  radius  u  lies  entirely  inside 

u>.  Multiply  (3.13)  by 


4  (u,  -  n)  (u.  >  n) 
h  h 


<u.  <  n) 

h 


and  integrate,  but  in  case  h  =  n  restrict  n  >  0.  Integrating  by  parts  we  have 

3u  ,  - 

/  <J$.  .  ^  +  6.  )  rf-  [(  (u.  -  M)JdC  =  0  .  (3.19) 

‘ -  1]  dC .  1  at.  h 


wn{u>n} 


Since  there  is  no  surface  integral  in  (3.19),  we  only  need  to  check  the  case  i  =  n. 
When  h  =  n,  on  f  by  (3.18)  we  have  u^  -  n  <  0.  When  h  ^  n,  j  =  n  from  (3.18)  we  have 


3Vr  3Vr 


—  =0.  When  h  n,  j  ?  n,  from  (3.11),  (3.14)  we  have  8.1  =  0.  And 


the  b,  defined  by  (3.15)  has  the  property  b  )_  »  0  when  h  +  n,  the  reason  is,  from 
i  n  [  p 


(3.11)  we  have  o  I  =  0  (j  f  n)  ,  and  this  induces  rrr- 
nllf  3f-r 


=0  in  the  case  j  K  n, 


h  f  n,  combine  with  (3.18)  we  have  b  J  *  0. 

n|r 

From  (3.19)  by  applying  Schwarz  inequality  and  (3.17)  we  have 


k,u-8u 


grad  u,  )  dt  4  L„  mes(8  n  A  )( - -  max  [u. 

^  9  n'u  (Ou)2  te5nA  ^ 


[u.  (t>  -  Ml  +  1}  .(3.20) 


/VV - Cn-1' 

; . 


-  C  )  (h  t  n) 
n 


(0  )  -  1  - 

in  the  part  of  to  outside  S2 ,  then  apply  (3.18)  we  have  u  e  C  U)  W  .,(»)•  from 


(3.20),  (3.21)  we  have 


/  (grad  u.  )2d£  4  2L_  mes  A  { — — -  max  [u  (?)  -  n]2  +  1} 

n  y  Hip  */\  »  n 

A  „  '  0u  CeA 

n,u-9u  n,u 


(restrict  n  >  0  when  h  -  n). 

By  the  same  process  we  can  obtain  two  formulas  similar  to  (3.4),  (3.7). 

From  the  result  of  [3]  we  have,  u^(h  =  1,...,n)  satisfy  a  Holder  condition  in  uj 

r2 


i.e.  there  exist  L10  and  t2(Q,  ft)  (0  <  t2  <  1)  such  that  for  any  x,  y  e  ft  n  w  we  have 


|u(x)  -  u(y) |  4  L1q |x-y| 


(3.22) 


§4.  Estimate  of  solution  in  the  neighborhood  of  «• 


Assume  that  the  solution  of  (1),  (8),  (9),  (11),  (2.17)  satisfies 

e  c(2) (ft)  n  c(1 * (ft  +  r>  . 


(4.1) 


When  (  £  E,  from  (2.5)  and  applying  the  Schauder  interior  estimate  twice  we  have 

•f  e  C(3+Tl(ft)  n  C(1,(ft  +  T )  . 


Let  O<041>O<S<1(A  is  a  positive  constant,  R,  R  are  large  positive 
numbers.  Let  r 


-  C(x,  R,  R)  *■  l  R  (1+0) 


.A  |x|-R 
Or 


R(  1+0  )  <  | x |  4  R(1-9  ) 


R  4  | x )  4  R( 1+9 ) 


RX(1-0  )X  R~ LX1  R(1-0  )  4  |x|  4  R  . 

V 


(4.2) 


Using  this  C  we  can  obtain  (3.3)  also.  When  X  <  1  -  — ,  from  (2.17)  and  (4.2)  we  see 
that  (3.3)  is  also  true  when  C  *  C(x,  R,  °°). 


Denote  the  inversion  point  of  x  about  the  unit  sphere  by  X,  then  (3.3)  becomes 

2.  ,4-2n,A  _ .2 


/  (;2|x|4*2n(Axu)2dx  4  l2  /  u':|X|‘l'^n(AxC)‘:dX  . 

Take  X«2-n,9*1,R*«’  and  let  —  «  U,  from  (4.2)  and  (4.3)  we  have 

/  (grad  u)2dX  4  L  un  2 
|X|4*  X  11 


Similarly  we  have 


2  1  2 
( urad  u)  dX  4  L  mes  A  - -  max  (u(X)  -  n] 

A  ,,  ,  -1  n,M  (0U)  XtA 

n ,u ( l+o  )  n  ,  u 


(4.3) 


(4.4) 


(4.5) 
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B 


(grad^  u)  dX  <  mes 


n#  u( i+S) 


- -  max  [u(X)  -  T)]4 

(9 u)  XeB 

n ,  u 


(4.6) 


where  A  ,  B  are  the  sets  of  u(X)  >  n  in  I X 1  <  u  respectively. 

h  #  U  0#  U  J 

It  is  easy  to  prove  that  u(X)  t  in  a  neighborhood  of  the  origin  (which  includes 

the  origin).  From  (2.17),  (4.4),  (4.5),  (4.6)  and  applying  the  result  of  [3],  we  see  that 
constants  a  -  a(Q,  D  >  0  and  L14  exist  such  that  in  the  neighborhood  of  the  origin  the 
following  inequality  is  valid: 

|u(X)  "  u|x=0l  4  L14|X|C1  (4.7) 

where  ulx=Q  =  ulx=<„  =  u<“>  is  defined  from  (11).  Let  u(x)  -  u(~)  =  v(x),  then  (4.7) 
becomes 


I v ( x ) |  <  L14lx| 


(4.8) 


Apply  a  linear  transformation  x  =  x(y)  such  that  the  matrix  (a  )  becomes  an 

lj  x=» 


unit  matrix.  Then  (3.1)  becomes 


3  ,  3v  .  „ 

3vi  ^  3yj 


where  c.  .(»)  -  5.  .  .  Combining  with  (4.8)  we  have 
1 J  i  1 


\o.Ay)  -  «..|  <  L15lyl'a  . 


2 


Multiply  (4.9)  by  ;  'v  and  integrate  to  obtain 

/  52(grad  v)2dv  <  L  f  v2 ( grad  ;)2dy 


(4.9) 


(4. 1(1) 


(4.11) 


by  us'.iq  the  same  method  as  (3.3)  was  obtained,  where  t,{  y)  =  <;(y,  R,  R)  is  determined 
from  (4.2)  by  changing  x  to  y.  Take  8  satisfying 

0  <  8  <  j  -  1,  B  a  .  (4.12) 

In  (4.2)  take  \=1-  —  +3,  R  =  «>,  9=1,  substituting  in  (4.11)  and  applying  (4.8)  we 


f(r)  =  f  |y|2\grad  vl^dy  4  L  r 

ly  I  >r 

Estimate  bv  extending  the  method  of  [41.  Denote  the  unit  spherical  surface  with 


2(  B-a) 


(4.13) 


anter  origin  by  let  r  <  R,  and  denote  1  y  1  =  y  we  have 
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R  2 

f(r)  -  /  |y|2X(grad  v)2dy  >  /  dw  /  (|^)  pn  1+2Xdp 


r< I y I <R 


(4.14) 


R  R  2 

>  f  [/  ^  dy/ f  u  n+1  2Xdy]  du  >  (n-2+2X)rn  2+2X  J  (v(B)-v(r) ]2du 


where  in  the  above  expression  Schwarz's  inequality  was  used.  Let  R  ♦  »  in  (4.14). 


Applying  (4.8),  (4.13)  we  have 


f(r)  >  (n-2+2X)rn  "'f21  f  v2(r)dw  -  (n-2+2X)  /  I  y  1 2  ^  1v2dS 


lyl-r 


Combining  (4.12),  (4.13),  (4.15)  we  have 

r 

I  y  I  >r 


,  ,  .28-2  2  .  _  T  2(8-0) 
/  I y I  v  dy  <  Llgr 


Let 


g  ( r )  «  /  I  y  I  2  Xc  .  .  dy 

R  r  <  |  y  |  <R  13  3yi  3yj 


Integrating  by  parts  and  applying  (4.9)  we  have 


g  (r)  =  (J  -  !  )|y|2Xvc  —  cos(N,  y  )dS  - 

R  I y l=R  lyl-r  3yi  3 


(4.15) 


(4.16) 


(4.17) 


(4.18) 


!  1  y I  2 X  'vc  ~  cos(N,  y  )dy 

r  <  |  y  |  4R  *3  3yj  3 


where  N  is  the  normal  directed  toward  «.  When  R  =*=  «>  from  (4.10),  (4.13)  we  see  tht 
the  left  hand  side  of  (4.18)  has  a  limit,  from  (4.13),  (4.16)  applying  Schwarz  inequality 
we  prove  easily  that  the  last  term  in  the  right  hand  side  of  (4.18)  has  a  limit  when 
R  ♦  ->,  hence  when  R  ♦  •» 

f  ly|2Xvc  cos(N,  y.  )dS  (4.19) 

|  y  |  =19  ^  ^i  3 

has  a  limit.  If  this  limit  value  is  a  constant  p  ¥  0,  then  applying  (4.10),  (4.13), 


(4.15)  when  u  large  enough  we  have 

p2  <  L  (f  ( y  I  2  x  I  v  |  I  grad  v|ds)2  <  L  u  f  ly|2>‘  1  v2dS 


lyl=u 


lyl  =u 


Integrating  we  get 


f  |y|2X(grad  v)2dS 

I  vl=u 

<  l-t  liflulf'llilflull/hi  -  2  +  2X)  . 


p2(n-2-2X)1nu/L1(1|R  «.  -  ^  f2(y>l!?  (R0  <  R) 

n  o 


(4.20) 
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(4.13)  contradicts  (4.20)  when  R  +  »,  hence  the  limit  of  (4.19)  must  be  zero.  Let 

R  +  "  in  (4.18)  and  denoting  qjr)  by  g(r)  we  have 

g(r)  *  -/  ly|2Xvc  |~7-  cos(N,  y  )dS  -  2X  J  ly|2X_1  v  dy 
1 Y I =r  1  yl  3  ly | >r  j 


-2X  / 


ly|2X  v  ( c  ,  -  5  .)  -j— —  cos  (N ,  y.)dy  =  1+1+1 

|y|>r  15  1]  3  12  3 


Applying  (4.10),  (4.13),  (4.17)  we  have 


(4.21) 


g(r)  =  (1  +  0 ( r  “)lf(r> 


(4.22) 


Applying  Schwarz  inequality  to  I ,  and  applying  (4.10),  (4.13),  (4.14)  we  have 

x2  <  r  /  ly|2X  1v2ds  /  lyl  2X[c  cos(N,  y  )  ]  2ds  4 

I  yl=r  |yl=r  ij  3yi  j 


SSk/ ..  lv!2Xt-|+  lvra0,|gradv,,2dS  4 

1  y  I  =r 


[1  +  0 ( r”a) ] rf ( r )/(n-2+2X)  j  |y|2X(grad  v)zds 

lyl=r 


-11+  0( r  j] [rf(r)f,(r))/(n-2+2X) 


(4.23) 


Integrate  by  parts 


T  _  r  2X+n-2.  ,  sv  ,  ,  2A+n-2  ,  2,  , 

I2  '  "2X  j  u  dp  j  v  —  da)  «  -X(p  J  v  dui) 


3v 


2  , 


3U 


U=Y 


+  X( 2X+n-2 )  /  u2X+n"3dg  J  v2du  »  <1  -  |  +  S)  f  |y| 


2X-1 


«  I y I =r 

v2dS  +6(1  -  ^  +  6)  /  lyl 2X_2  v2dy  4  0 

lyl-r 

the  validity  of  the  above  formula  has  used  (12)  and  the  relation 


(4.24) 


From  (4.15)  we  have 


,  2  X+n-2  ,  2 , 

(  g  f  v  dw)  -  0  • 

‘  g“» 

0) 


lyl  4r 


.  . 2X-2-2g  2  .  ,  1 

1 Y I  v  dy  < 


n-2+2X 


/  u  1  2ady  /  |y|2A(grad  v)2dy 


2X, 


ly  l=r 


<L20r-2af(r, 


(4.25) 


Prom  (4.10),  (4.13),  (4.25)  and  applying  Schwarz  inequality  to  estimate  13  we  have 


|xg'  4  Si  / 


j  ly|2X  1  a|v|  I  grad  v|dy  4  L  r  "“f(r) 

ly I >r  22 


(4.26) 


Substituting  (4.22),  (4.23),  (4.24),  (4.26)  into  (4.21)  we  have 


[1  +  0  ( r_a)  1  f  (  r )  4  ((1  +  0(r'a)1  (-rf(r)f*  (r)  ) /(n-2+?\)}/2 
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Integrating  we  get 


f(r)  <  [1  +  0(r  “)) <-rf' (r)1/(n-2+2A) 


r  n-2+2A  P(r'“) 


f(r)  <  f  ( r  )  e  0  <  L,,r2"n‘2X  (r  >  r  (Q,  D  ) 

0  '  23  o 


,  ,  ,2-n+28,  „  >2,  .  -2B 

)  |y|  (grad  v)  dy  4  L  r 


lyl>r 

Substituting  (4.27)  into  (4.15)  we  have 


Take  y  satisfying 


/  lyl‘n+28v2dS  4  L  r” 
lyl-r 


0  <  Y  <  j  -  Y  <  20  . 

In  (4.2)  taking  Xm1-~+Y#R“00»  0  «  1,  substituting  into  (4.1)  and  applying  (4.28) 


/  I  y  I  2  n+2'((qrad  v)2dy  4  (;2(grad  v)2dy  4  t_.  f  I  y  I  n+2lfv2dy  4 


T  ,  -1-4B+2y  -40+2y 

L26  U  <  L27R 

Note  that  (4.29)  improves  (4.13). 

Repeating  the  above  process  we  prove  that 

f  I y! 2  n+2>(grad  v)2dy  4  I.  r  2y  . 

I  y  I  >r 

Continuing  in  this  way,  after  a  finite  number  of  steps  we  have 

/  (grad  v)2dy  4  L  r2  n  . 

ly  I  >r 

We  need  the  followinq  result  of  [21  extended  to  the  n  dimensional  case: 

/  (grad  v)2  dy  <  L  r  2*^_1 

I  Vl  >r 

Proof.  Let  S  be  a  closed  surface  which  lies  in  ft  and  contains  r  in  its 
interior.  Integrating  (1)  and  applying  (9)  we  have 


f  0  cos(N,  x .  )  dS  *=  f  OT”dS  =  n 
-  dx .  l  1  _  3N 


When  x.  changes  to  x  t  e,  the  chanqed  function  will  be  denoted  hv  adding 
n  h 

subscript  c.  Denote  bv  S  the  surface  5  translated  bv  (n,...,0,  -r ,  n,...,h). 


from  (4.32)  we  have 


/  (p  -r — )  cos (N/  x.)dS  ■  f  p  — —  cos(N,  x  )dS  *  0 

I  3Xi  E  1  S  3xi  1  E 

e 

Differentiating  with  respect  to  the  parameter  t  and  let  c  »  0  we  have 


I 


f  »  -5-Se- 


aij  it cos(N'  Vds  =  L  aij  if-  cos(N'  Vds 


s  J  i 


's  3xi 


f  fp  I  „  cos(N,  x  )dS  »  0. 

3e  3x  'e| e=0  i 


After  linear  transformation  (S  changes  to  S)  we  have 

/  c  cos(N,  y  !dS  -  0  .  (4.33) 

S  3  dyi  3 

From  (4.30)  we  see  that  we  can  take  X  *  0  in  (4.21),  and  it  become  (by  using  (4.33)) 
g(r)  =  -f  (v  -  v)c.  .  cos(N,  y.)dS  (4.34) 

'  .  .  It  4  t 


v  =  J  vdS/mes  S 

I  y  I  =r 


From  wirtinaer's  inequality 


f  (v  -  Z)2dS  i  ~  /  [(grad  v)2  -  (|~)2ldS 

ly|-r  lyl-r  ™ 

and  (4. in),  (4.34)  we  have 
/n- 1  , 

9  r  <  2r  '  .  ,  ( v-v)  dS  +  /  [c .  .  cos (N,  v  . )  ]  2dS  < 

1 vl=r  2/H^T  'ly|=r  n  3y.  j 


27h=t  ,  f(qrad  v)2  -  'ds  * 


Y7=f  /  [  ^r  n(  grad  v)2|ds  =  [1  +  0(r  “> )  f  (grad  v)2ds 

lyl  =  r  '  '  I  y  I  =r 

=  -r  [  1  *  0  ( r”a)  ]  / (  2/n-T )  g'  ( r )  . 


Inteqratinq  we  have 


.r0.2/n-1  0(r0  1 


g(  r  )  <  g(  r0  )  <  — )  e 


(r  >  r0(Q,  D)  . 


From  (4.10),  (4.17),  (4.30)  we  obtain  (4.31). 


Combining  (4.30),  (4.31)  we  obtain 


lyl>r 

Substituting  into  (4.15)  we  have 


(grad  v)2dy  4  L^r"  4  (6  >  0) 


lyl«r 


2  .  n-3-26 

v  ds  <  L33r 


Hence  we  have 


,  i  1 2-n+6  2,  ,  2-n+  6 .  , 

J  I y |  vdy-Jr  dr  J 


I  y  I  >R  R 

Reverting  to  the  variable  x  we  have 


lyl«r 


v2dS  4  L,  R*{ 
34 


/  I y 1 2  n+^  V2dx  4  L  _  R  ^  .  (4.36) 

| x | >R  35 

Let  \Ji  «  *-Ux  where  U  is  defined  by  (11)  above,  then  from  (4.36)  we  have 

/  ly  I  2  n+^(grad  i|;)2dx  4  nL  -R  ^ 

|x|>R  35 

Let  X  be  the  inversion  point  of  x  about  the  unit  sphere  and  let  —  »  u»  the  above 
formula  becomes 

/  |X|2_n'i(gradv  *)2dX  4  nL,_uS  •  (4.37) 

0  |X|<U 

Fixing  a  point  x  e  0,  without  loss  of  generality  we  can  assume  temporarily  that 
«o  0  x 

*(x°)  *  0,  hence  \Mx)  “  -q  x1  +  /  Q  vhdx^,  and  combining  with  (2.17)  we  have 

x 

I *■!*>!  <  L^lx] 


U(X)|  4  L-  I  X I 
Jw 


-1 


(4.38) 


From  (4.37),  (4.38)  we  have 


*  €  W2  in  I X |  4  u0(C,  D 


and  the  following  Sobolev  decomposition  formula  is  valid  [12]: 

lMx)  =  /  C<v)'Mx)dy  +  /  |y|1_nw  (x,  Y)  dy 

|y|<u0  i y i <u0  3Yh 

where  £  and  w.  are  known  bounded  functions. 


(4.39) 


From  (4.37)  we  have 
.  1-n , 


f  I  Y|  |  grad  *|dY  <  [/  |Y|0’ "dY  / 

|y|<u  I y I <u  ' l y l 


i  2-n- 6, 


t  2hv  1V2 


Hence  the  right  hand  side  of  (4.39)  ig  a  continuous  and  bounded  function  when  |x|  <  . 


Therefore  0(*)|  exists,  and 

x*® 


f  '!>(  x  )  (  <  L 


39 


(4.40) 
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Denote  t |i(x)  -  I still  by  ^(x).  Then  from  (3.3),  (3.22),  (4.8),  (4.40)  we  obtain 
that  the  solution  of  (1),  (8),  (9),  (11),  (2.16)  is  bounded  in  E. 

All  the  requirements  of  Leray-Schauder  degree  theory  are  satisfied,  hence  theorem  2  is 
proved. 


§5.  Some  properties  of  subsonic  flow. 

Theorem  3.  The  solution  of  problem  (1),  (8),  (9),  (10)  is  unique  in  the  subsonic 
range. 

Proof.  From  (2)  we  obtain  that,  the  difference  y  -  v>^2*  -  of  two  solutions 


>  satisfies  the  following  equation 
(2) 


a..(u'-' . u(2),^|-+ 

11  1  n  JxJ*.  k  dx, 

ilk 


(5.1) 


,  ( 1  )  o\D  (2)  a'fi  . 

where  u,  *  —5 - .  u.  =  -  and 

i  dx.  i  dx 


b  =  [a  ( u 
k  1  i  j  1 


(2) 


(2)  (1) 

,Uk  '  Vi' 


( 1 )  ,  ,  (2)  (2) 

'Un  )_aij(ui  . Vi' 


(1) 


2  (  1  ) 

(1).,  ,,  (2)  (1),  3 

,Un  ,1/tuk  '  Uk  1 


(5.2) 


From  (8)  we  have 


(5.3) 


From  the  last  section  we  have:  ^2^  (after  subtraction  of  a  suitable  constant) 

satisfy  (2.6),  hence  we  have 

£(«)  «  0  .  (5.4) 

From  (5.1),  (5.2)  by  applying  the  maximum  principle  [5]  we  have,  when  7  is  not  a 

constant,  then  it  cannot  taken  a  positive  maximum  value  or  a  negative  minimum  value  in 

+  r,  so  combining  with  (5.4)  we  have  7  “  0.  This  proves  the  theorem. 

00  00 

Theorem  4>  There  exists  a  positive  constant  q^#  such  that  when  0  <  U  <  q^,  the 
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solution  of  problem  (1),  (8),  (9),  (10)  exists  in  the  subsonic  range  (this  solution  is 


unique  by  theorem  3).  Moreover,  the  function  Q  defined  by  (12)  satisfy 
Q (U )  e  C  [0,  q^) ,  and 

lig  Q(U>  -  q  .  (5.5) 

U+qc-0 

The  proof  is  divided  into  several  points  as  follows: 

1)  Let  }  be  any  sequence  of  solutions  of  (1),  (8),  (9),  (11),  (12),  (13)  with 

n 

Q,  U  substituted  by  Qn,  On#  If  sup  Qn  <  qc,  then  from  §2,  iii  we  have  that  11*9  I  is 

uniformly  bounded  in  E,  is  compact  by  theorem  1,  hence  a  limit  function  *p  exists, 

it  is  easy  to  prove  that  $  is  a  solution  of  (1),  (8),  (9),  (11),  (12),  (13)  also  and  its 

corresponding  Q  is  a  limit  point  of  the  sequence  {Q  } . 

n 

2)  Take  any  solution  f  of  (1),  (8),  (9),  (11),  (12),  (13),  we  define  U  «  U(Q,y') 

by  (11).  Fixed  Q(0  <  Q  <  q^)  and  denote  the  infimum  of  all  U(Q,V>)  by  V(Q),  i.e.  a 

solution  sequence  }  exists  such  that  U(Q,V  )  *  V(0).  By  1),  there  exists  a  limit 
n  n 

function  »  ^(Q)»  which  is  the  solution  of  (1),  (8),  (9),  (11),  (12),  (13)  (where  we 
substitute  for  U  by  V  in  ( 1 1 ) ) . 

3)  Let  0  <  <  q^,,  we  prove  that  the  interval  [0,  V(Q0 )  ]  is  covered  by  the  set 

of  all  (U,  v>)  (0  <  Q  <  Qq)  . 

If  this  is  false,  i.e.  constant  Ug  exists  such  that  0  <  Ug  <  V(Qg)  and  the  problem 
(1),  (8),  (9),  (10)  (where  in  (10)  we  substitute  for  U  by  Ug),  has  no  solution  in 

0  <  Q  <  Qg. 

If  there  is  a  positive  number  sequence  e  with  e  ♦  0  ( r»  ♦  00 )  such  that  (1),  (8), 

n  n 

(9),  (10),  (where  in  (10)  we  substitute  for  U  by  Un  -  e  ) ,  has  a  sequence  of 

u  n 

solution  yfi n  in  0  <  Q  <  then  by  1),  the  limit  function  is  the  solution  of  (1), 

(8)/  (9),  (10),  where  in  (10)  we  substitute  for  U  by  Uq,  in  0  <  Q  <  Qq,  and  this 
contradicts  the  above  hypothesis# 

Hence  a  constant  U^  exists,  satisfying  0  <  <  Uq,  such  that  the  problem  (1), 

(8),  (9),  (10)  has  no  solution  in  0  <  Q  «  QQ  when  U^  <  U  *  U^# 

l'\ 


Now  apply  Leray-Schauder  degree  theory*  to  the  set  -<•  <  U  <  V2  (UQ  +  U  )  in  the 
space  E  (for  the  relation  of  the  element  •f  e  E  and  U  see  (2.6)).  We  obtain  a 
solution  to  problem  (1),  (8),  (9),  (11),  (12)  for  any  given  0  <  Q  <  Qq  and  this  solution 
possesses  the  following  property:  V  defined  from  (9),  (10)  satisfies  0  <  U  <  UQ.  Hence 

inf  V(Qq)  <  V(Q0)  . 

This  contradicts  the  definition  of  V(Qg),  which  proves  3. 

4)  The  function  V(Q)  is  strictly  monotone  increasing  in  the  interval  (0,  q^. 

Otherwise  there  are  0  <  Q,  <  Q.  <  q  such  that  V(Q  )  >  V(Q, ).  From  3)  we  have,  when 

12c  12 

0  <  Q  <  Q1  the  set  of  U(Q,^)  covers  the  Interval  [0,  VfQ^)],  hence  we  have 
Q0  £  [0,  and  a  such  that  U(Qn,  *fg)  “  V(Q2),  contrary  to  theorem  3. 

V(o)  is  strictly  monotone  increasing  and  has  supremum  qc,  hence  V(qc  -  0)  exists. 

5 )  Let 

q"  =  V(q  -  0)  . 

c  c 

From  3)  we  have,  the  set  of  U(Q)  (0  <  Q  <  q  )  covers  the  interval  0  <  U  <  q*.  In  other 

c  c 

words,  the  solution  of  problem  (1),  (8),  (9),  (10)  in  the  subsonic  range  exists  when 

0  <  U  <  Q°°» 
c 

The  above  solution  is  also  unique  by  theorem  3,  hence  the  function 

Q(U),  0  <  0  <  q"  , 

is  determined  uniquely. 

For  fixed  U0  in  (0,  q^) ,  it  lies  in  the  number  set  U(Q(Uq),1)0  hence  we  have 

V(Q( UQ)  )  <  UQ  .  (5.6) 


*  . 

From  §2,  i i i 0  but  substitute  Q  by  positive  constant  K  exists  such  that  the 

solution  of  (2.13)  satisfying  j  |  v  ||  K.  Take  F  be  the  part  of  sphere  ||tf||  <K+1 
satisfying 

-  w  <  U  <  V2  (Un  +  Ut  )  =  Ux  t  +  *  ♦  ,  D  V t  ( oc )  =  0  ) 

then  (2.13)  has  no  solution  in  the  neighborhood  of  the  boundary  of  F.  From  §2,  ii  we 
have,  the  lecrree  of  solution  is  1  in  F  when  *  =  n,  combine  with  ^2,  {  ^nd  applying 
Leray-Schauder  legree  theory  we  obtain  that,  when  0  ^  1 ,  (1)#  (B),  (9),  (11),  (12) 

In  F  has  at  least  one  solution. 
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6) 

From  the  definition 

of  q°°  there 
c 

exists 

«n 

♦  q  -  0,  such 
c 

that 

v ( Q_  >  *  q“ 

n  c 

If 

u  satisfies  0  <  u 

<  V(o  )  and 
'  n 

Q(U> 

> 

V 

Combining  with 

(5.6) 

we  have 

U  >  V(Qn 

) ,  this  contradicts 

U  <  V(Q  ). 
n 

Hence 

when 

0  <  U  <  V(0  ). 

n 

we  have 

0  <  Q(U)  <  Q  <  q  .  (5.7) 

n  c 

7)  Let  {U  }  satisfy  0  <  U  <  V(Q  ) ,  0  ♦  U  ,  Q(U  )  ♦  Q  .  From  (5.7)  by  applying 

in  m  n  m  u  m  o 

1)  we  obtain  that  the  limit  function  f  is  the  solution  of  (1),  (8),  (9),  (10),  where  U, 
Q  are  given  by  Ug,  Qg  in  (10),  (12)  respectively.  Applying  theorem  3  we  have 
Q0  “  Q(Ug),  in  other  words  Q(U)  c  C[0,  VtQ^)].  Letting  n  ♦  «  we  have 

fi(U)  E  C[0,  q“)  . 

8)  From  the  above  points  we  see  easily  that  lim  Q(U)  »  q  .  then  the  functional 

U+q“-0  C 

c 

value  of  the  continuous  function  Q(U)  oscillates  finitely  in  the  interval  (Qg,qc)»  and 

every  point  in  (Q0,qcl  is  the  limit  point  of  Q ( U )  when  U  ♦  q  -  0.  Hence  we  have 

U  *  q  ,  U*  ♦  q  such  that  Q(U  )  ♦  Q_ ,  Q(U')  ♦  V?  (Q„  +  q  ).  By  1),  the  limit  functions 
mcmc  mOm  uc 

f ,  if'  of  corresponding  solution  sequences  fa  },  fa ' }  are  solutions  of  (1),  (8),  (9), 

m  in 

(11),  (12),  (13),  where  Q  is  given  by  QQ  and  Vj  (Q0  +  qc>  respectively.  Hence 

„ 

O(Q0)  -  qc  >  0(-Y~) 

or  the  problem  (1),  (8),  (9),  (10)  (substitute  U  by  q^)  has  at  least  two  solutions  in 
the  subsonic  range,  which  contradicts  theorem  3,  hence  (5.5)  is  true. 

The  theorem  is  proved  completely. 

Condition  (11)  restricts  the  incoming  flow  is  in  the  positive  x1  direction.  Now 

remove  this  restriction  by  only  assuming 

00  00  00 

u=(u.,...,u)  (5.8) 

I  n 

is  given.  Because  (1),  (8)  is  invariant  under  axis  rotation,  hence  by  theorems  3,  4  we 
have,  the  solution  of  (1),  (9),  (9),  (5.8)  exists  and  is  unique.  Let 

> p  =  u^  +4)  (5.9) 

and  regard  as  a  solution  of  the  linear  nroblem,  by  theorem  1  we  have 


,*,n-2.2+t  4  V 


(5.10) 
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Now  prove  that  the  solution  of  (1),  (8),  (9),  (5.8)  depends  continuously  on  u  by  the 
following  sense. 

Theorem  5.  There  exists  a  positive  constant  L4n,  such  that  for  any  two  vectors  u°°' ' , 

oo  2 

u  '  in  the  class 


I  u  |  <  V(-* —  u,,  uj 

lu  |  lu  | 

(the  definition  of  function  V  see  theorem  4)  the  corresponding  solutions  •?'  ,  0f 


(1),  (8),  (9),  (5.8)  satisfy 


tl<2)  ,(1). 
!  4  -41  I 


06  2  CO  1 

„  „  4L..U  -u 

n-2, 2+t  40 


(5.11) 


where  the  relation  of  and  4  is  given  by  (5.9). 


Proof.  Let  *  i|/2'  -  i/1*,  from  (2)  we  obtain  that  ?  satisfies  (5.1),  (5.2). 


Prom  (8)  we  have 


•^Ir  =■ -(u"'2  -  U*'1)  cos(N,  Xi)!r  . 

Applying  (6),  (5.10)  it  is  easy  to  show  that  the  b^  defined  from  (5.2)  satisfy  conditions 
(1.7),  (1.55),  (1.59),  and  from  (5.10)  we  get  (5.4).  Hence  from  theorem  1  we  have  (5.11). 
This  completes  the  proof. 

Theorem  6.  For  every  non  zero  subsonic  flow  around  a  given  profile,  cannot  be  taken 

on  1 1  or  at  ■».  In  other  words,  qmax  can  only  be  taken  on  P.  And 

<W<L41|U"'  ’  <5-12) 

Proof.  If  >  0  is  taken  by  P  c  Cl,  then  after  a  rotation  of  axis  we  arrive  at 

grad  u|p  =  ^max'  °*«**»°)  and  (1),  (8)  remain  unchanged,  hence  u^  takes  the  maximum 

value  at  P.  And  u1  satisfies  the  elliptic  equation  (3.1),  hence  from  the  maximum 

princide  (51  we  have,  \iy  is  a  constant,  i.e.  =  amax,  hence  u^  -...=  un  =  0.  From 

(8)  we  have  a  =  0.  This  is  a  contradiction, 
max 

Prom  (5.10)  and  aoplying  the  result  of  (2J,  we  have  that  draax  cannot  be  taken  at 


M  2 

(5.12)  is  a  special  case  of  (5.11)  with  u  =0.  The  rr->of  of  the  theorem  is  thus 
complete . 

Theorem  7.  u)  exists  for  any  subsonic  flow  in  the  neiahborhood  of  «.  -»nd  (4.8)  is 

valid.  An  the  subsonic  flow  in  the  whole  space  is  uniform. 
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Proof.  For  subsonic  flow  in  the  neighborhood  of  «,  because  (4.4),  (4.5)  and  (4.6) 
are  valid,  hence  the  existence  of  u|  and  the  validity  of  (4.8)  follow  from  [3). 

X*® 

Concerning  subsonic  flow  in  the  whole  space,  u|  ^  exists  by  the  above  argument, 
hence  it  is  a  special  case  of  flow  around  profile,  i.e.  (4.32)  is  true  for  every  closed 
surface.  Hence  (5.10)  is  valid,  and  we  can  apply  theorem  6,  i.e.  except  q  is  a  constant, 
cannot  be  taken  by  any  finite  point  and  •>.  And  when  q  is  a  constant,  u  is  a 
constant  vector  by  theorem  6,  i.e.  the  flow  is  uniform. 

The  theorem  is  thus  proved. 
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